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ABSTRACT 


An  approach  to  prediction  of  long -period  Rayleigh  wave  signal 
persistence  is  taken  through  "mean-time -between-failure"  and  "waiting  time" 
stochastic  models.  The  estimate  applies  only  to  the  generalized  case  of  a 
network  of  recording  stations  and  wide  spatial  distribution  of  earthquake 
sources. 

Rayleigh  waves  from  Eurasian  earthquakes  recorded  on  Very 

Long  Period  Experiment  systems  were  analyzed  to  obtain  data  for  the  study. 

Successive  time  intervals  between  the  maximum  recorded  amplitude,  A  ,  and 

the  latest  point  of  recorded  signal  equal  to  A/ 2,  the  time  between  A/ 2  and 

Aj 4,  A/4  and  A/ 8,  and  etc.  were  measured  until  the  signal  amplitude  reached 

the  ambient  noise  level.  Resulting  "half -amplitude  decay  times"  (W  )  were 

n 

found  to  follow  a  gamma  probability  law.  The  distribution  of  observed  W  fits 
this  law  with  0.  90  confidence  or  greater  according  to  the  Chi- Squared  criterion. 
The  results  indicate  that  the  probability  distribution  of  half- amplitude  decay 
times  in  earthquake  Rayleigh  waves  can  be  readily  estimated  using  appropriate 
gamma  distribution  parameters. 
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SECTION  I 
INTRODUCTION 

Estimation  of  the  significance  of  seismic  signal  interference 
in  the  operations  of  a  network  of  seismic  recording  stations  is  of  importance 
both  toward  understanding  signal  detection  potentials  and  for  evaluating  tne 
necessity  of  development  and  use  of  "anti-interference "  analytical  tools.  Use 
of  the  seismograms  from  the  Very  Long  Period  Experiment  (VLPE)  stations 
provides  the  opportunity  for  study  of  generalized  multiple  path  signal  trends 
in  a  network  sense.  Current  seismological  trends  are  toward  Generalization 
of  operations  to  the  point  where  world-wide  station  distributions  are  treated 
simply  as  instrument  arrays,  and  the  kind  of  approach  taken  here  is  both 
necessary  and  vital  to  successful  use  of  data  provided  by  the  recording  stations. 

While  most  who  are  familiar  with  seismograms  can  call  to  mind 
some  approximate  length  of  time  that  long -period  Rayleigh  waves  are  visible 
on  our  recordings,  virtually  no  numerical  data  on  the  subject  have  been  pub¬ 
lished,  and  no  time -amplitude  persistence  reports  are  known  to  us.  The  sub¬ 
ject  is  of  interest  upon  occasion,  particularly  if  we  wish  to  have  a  measure 
of  some  specific  signal  at  a  time  when  another  signal  occupies  the  seismogram. 
The  length  of  time  that  some  level  of  amplitude  of  the  "unwanted"  signal  per¬ 
sists  is  also  of  some  interest  if  we  wish  to  forecast  the  possibility  that  the  sig¬ 
nal  interference  is  at  a  level  where  the  required  measurement  will  be  obscured. 
We  hope  to  provide  a  method  for  estimation  of  the  time -amplitude  characteris¬ 
tics  of  long-period  Rayleigh  wave  signals  with  this  study,  using  signals  record¬ 
ed  by  the  VLPE  stations  from  Eurasian  earthquakes  as  a  data  base. 


SECTION  II 
THEORY 


Stochastic  models  which  describe  emission  of  radioactive 
particles,  the  flow  of  telephone  or  automobile  traffic,  mean  time  between 
component  failure,  and  other  random  occurrences  of  events  in  time  are  wide¬ 
ly  used  to  provide  estimates  of  event  occurrence  probabilities.  The  "event" 
selected  here  for  analysis  is  the  occurrence  of  signal  amplitude  decay  to  one- 
half  its  maximum,  the  decay  of  that  half-amplitude  to  one-half  of  its  value, 
and  so  on.  The  time  required  for  this  event  to  occur  is  the  factor  we  wish  to 
describe.  We  therefore  need  to  find  a  model  which  will  fit  the  observed  data 
and  be  able  to  assign  numerical  values  to  the  model  parameters.  Such  models 
usually  fall  in  a  general  class  termed  "inter -arrival"  or  "waiting  time"  sto¬ 
chastic  processes  (Parzen,  1962.  1967;  Hogg  and  Craig,  1968). 

One  basic  model  for  waiting  times  is  that  used  for  events  of  a 
Poisson  type.  The  model  parameter  is  the  rate  of  occurrence  of  innependent 
events  within  successive  time  intervals  of  length  t  .  If  Nk  >s  the  number  of 
times  that  k  events  are  observed  within  the  interval  t  .  the  probability  P(k) 
that  exactly  k  events  will  be  observed  in  t  is  found  to  be: 


P(k)  = 


e-M(At)k 

k! 


where  A  is  the  reciprocal  of  the  mean  numbe.  of  events  observed  in  t  (i.e., 

A=  1/m  if  m  is  the  mean  of  observed  times)  or  the  event/unit  time  rate  of 

occurrence  (  =  E  N  (k)/t).  The  model  may  be  used  to  describe  the  expected 

time  for  an  event  to  occur  as  well  as  the  probability  of  an  event  occuring  within 
some  discrete  time  interval  in  a  succession  of  intervals.  The  probabi'ity  P(k) 
that  at  least  k  events  will  be  observed  in  a  time  interval  is 


P(k)  =  P(n)  +  P{n+1 )  +  P(n+ 2)  +.  .  .  +  P(n+m) 


(2) 


n+m 

■  s 


n- k 


-At  .  .  .k 
e  (At) 

k! 


n  =  k,  k  H ,  k+2 ,  .  .  . 


and  the  probability  P  [k]  that  no  more  than  k  will  be  observed  in  the  in 
terval  is 


P  [k]  =  1  -  P(k)  O) 

For  this  study  we  know  that  an  event  has  occurred,  i.  e.  ,  that 
half-amplitude  was  reached  in  our  observations  simply  because  we  have  a 
measurement  (or  several  measurements).  The  Poisson  process  lead  to  use 
of  the  exponential  probability  law  to  describe  our  data,  or  more  generally, 
gamma  probability  functions  related  to  the  Poisson  process.  If  the  observed 
half-amplitude  decay  times,  which  we  will  now  call  W^,  follow  a  Poisson  pro¬ 
cess  of  occurrence  at  a  mean  rate  of  A  events/unit  time,  then  the  occurrence 
of  the  rtH  event  (r  =  1,  2,  3,  .  .  .  )  obeys  a  gamma  probability  law  with  para¬ 
meters  r  and  A  .  This  law  has  a  probability  density  function  f(t)  of: 

f(t)  =  jr$ TT  ta0  (4) 

=  o  t  <  o 

The  waiting  time  (or  decay  time)  to  the  first  event  (from  peak  amplitude  A  to 
A/2,  or  A/2  to  A/4, .  .  .  )  then  follows  the  gamma  probability  lav'  with  para¬ 
meter  r  =  1  and: 

f(t>  =  Ae‘At  (5) 

For  t  >0,  the  probability  that  the  waiting  time  to  half-amplitude  will  be  less 
than  or  equal  to  t  (or  that  the  interval  t  will  include  at  leas  one  half¬ 
amplitude)  is  the  integral 


II- 2 


p 


(6) 


/ 


and  the  probability  that  it  will  be  greater  than  t  (or  that  the  interval  t  will 
hold  no  hdlf-amplitude)  is 


For  the  problem  of  describing  successive  half  amplitude  times  after  peak 
amplitude  A  ,  one  easily  evaluated  extended  form  of  (7)  is  given  by  Parzen 
(1967): 


P 


-At 

r  =0,  1, 2, 


where  the  probability  distribution  of  occurrence  times  is  described  for  the 
first  half -amplitude  A/2  at  r  =  0,  for  the  second  half- amplitude  A/4  at  r  =  1  , 
for  the  third  half-amplitude  A/ 8  at  r  =2,  and  so  on. 

In  the  following  sections  the  distribution  of  all  half- amplitude 
times  observed  -  i.  e.  ,  the  time  between  peak  amplitude  A  and  A/2,  between 
A/2  and  A/4,  etc.  -  will  be  tested  as  a  Poisson  distributed  random  variable. 

We  will  show  that  goodness  of  fit  tests  for  observed  and  predicted  distri¬ 
butions  have  high  probability  results.  These  results  indicate  that  the  are 
independent  exponentially  distributed  random  variables  of  a  Poisson  type 
simply  because  the  Poisson  process  is  characterized  by  this  fact  (see  Parzen, 
1962).  The  demonstration  that  this  type  of  distribution  is  satisfied  will  allow  a 
prediction  of  the  distribution  of  times  between  A  and  successive  amplitude  in¬ 
crements,  varying  only  the  parameter  r  with  A  held  constant  (Equation  8). 

A  model  such  as  equation  8,  if  fitted  successfully,  would  allow  prediction  of 
the  expected  time  of  signal  persistence  on  the  recordings  as  well  as  the  am¬ 
plitudes  expected  during  the  continuation  of  signal  and  coda. 
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SECTION  111 

DATA  AND  PROCEEuRES 


A.  SOURCE  OF  DATA 

Data  for  the  study  were  recorded  at  VLPE  station  listed  in 
Table  111-1  (also  shown  in  Figure  111-1)  during  January  through  March  of  1972. 
Paper  record  seismograms,  analog  conversions  of  digital  magnetic  tape  re¬ 
cordings,  and  bandpassed  digital  conversions  to  analog  seismograms  were 
used  to  obtain  measurements  of  signal  amplitude  and  time  data.  Typical  re- 
sponse  characteristics  of  the  recording  channels  are  shown  in  Figure  111-2. 
Bandpassing  of  the  digital  recordings  was  applied  in  the  frequency  domain, 
retaining  data  in  the  19.  4  to  42.  5  second  period  range.  Only  vertical  com¬ 
ponent  data  channels  were  analysed,  and  only  those  Rayleigh  wave  signals 
which  appeared  free  of  interference  from  other  earthquakes  were  used  since 
we  were  attempting  to  isolate  the  basic  decay  patterns  of  the  waveform.  Pre¬ 
liminary  Determination  of  Epicenter  (PDE)  lists  were  used  for  hypocenter 
parameters.  The  general  location  of  sources  studied  are  also  shown  in  Fig- 

ure  111-1. 

D.  ANALYSIS  PROCEDURES 

Long-period  Rayleigh  wave  signals  from  these  sources  were 
analyzed  for  "half-amplitude  decay  times"  by  the  following  steps: 

1.  The  time  T  between  the  maximum  recorded  Rayleigh  wave 

amplitude  A  (peak  to  peak,  regardless  of  wave  period,  and  the  latest  recorded 
signal  amplitude  equal  to  A/2  (also  peak  to  peak,  regardless  of  wave  period) 
was  measured  to  the  nearest  10  seconds  and  recorded  with  distance,  station, 
bodywave  magnitude,  and  recording  system  for  each  signal  analyzed. 


Ill-  1 


TABLE  III-  1 


VERY  LONG  PERIOD  EXPERIMENT  (VLPE) 
STATIONS  AND  LOCATIONS 


Station 

Designator 

Latitude 

Longitude 

Charters  Towers,  Australia 

CT  A 

20.  IS 

146.  3E 

Chiang  Mai,  Thailand 

CHG 

18.  8N 

99.  0E 

Fairbanks,  Alaska 

FBK 

64.  9N 

148.  0W 

Toledo,  Spain 

TOL 

39.  9N 

4.  0W 

Eilat,  Israel 

EIL 

29.  3N 

34.  5E 

Kongsberg,  Norway 

KON 

59.  6N 

9.  6E 

Ogdensburg,  New  Jersey 

OGD 

41.  IN 

74.  6  W 

Kipapa,  Hawaii 
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Long  Period  Experiment  (VLPE)  Station  Locations  and  Seismic  Source  Areas  for 
Amplitude  Rayleigh  Wave  Decay  Signal  Analysis. 


PVT  ATIVE  RESPONSE 


FIGURE  III -2 

Typical  VLPE  Station  Vertical  Channel  Response 
Characteristics 
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Times  for  successive  half- amplitudes  (T^  T  ,  etc.)  for  A/4, 
A/8,  etc.  ,  respectively,  were  also  measured  to  within  10  seconds  to  the  place 
on  the  recording  where  the  signal  waveform  or  definite  signal-related  coda 
could  not  be  distinguished  from  the  ambient  seismic  noise.  The  approximate 
time  of  signal  lors  in  the  noise  was  also  recorded  if  the  half-amplitude  signal 
step  could  not  be  obtained  at  that  point. 


Several  examples  of  the  measurements  are  shown  in  Figures 
II1-3,  and  111-4  for  a  signal  source  in  Iran  as  recorded  by  stations  in  Norway, 
Spain,  New  Jersey,  and  Alaska. 


C.  DISCUSSION 

The  first  time  interval  T]  ,  the  difference  between  T  and 
t2,  T2  and  T3,  and  so  on,  represent  the  "half- amplitude  decay  times"  (  or 
waiting  times  to  half-amplitude)  that  we  intend  to  model.  These  times 

are  a  measure  of  Rayleigh  wave  persistence  viewed  through  the  recording 
channels  which  may  have  characteristics  suitable  for  modeling. 

The  problem  of  trying  to  characterize  half-amplitude  decay  is 
readily  seen  in  Figures  III-5,  III-6,  and  I1I-7  where  the  observed  T  ,  T  ,  T 
and  so  on,  are  plotted  versus  their  respective  amplitude  increments1.  Obser¬ 
vation  of  time  to  noise,  Tr  ,  is  plotted  half-wa/  between  the  last  observable 
amplitude  increment  before  signal  loss  and  the  next  increment  in  the  figures. 
The  plotted  data  are  restricted  to  observations  of  less  than  2200  seconds  total 

between  peak  amplitude  and  observation  time  to  allow  some  detail  of  the  main 
bulk  of  data. 

While  not  obvious  in  the  display  of  data  in  Figures  I1I-5,  I1I-6, 
and  m-7,  there  are  no  consistent  linear  trends  in  the  signal  decays  related  to 
distance.  An  exponential  randomness  will  be  shown  later  which  apparently 
dominates  the  amplitude  decay  trends  at  all  distances.  A  weak  tendency  for 
more  rapid  A/2  decay  at  close  distance  is  evident,  and  most  of  the  A/8  and 
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FIGURE  III- 3.  Measurement  Procedure  for  Half- Amplitude  Decays  of  Long  Period  Rayleigh  Wave 
Signals.  Epicenter  35.  ON  51.  OE  on  February  3,  1972,  m,  =  5.  0  (C&C-S) 


FIGURE  HI-4.  Measurement  Procedure  for  Half- Amplitude  De~ays  of  Long  Period  Rayleigh  Wave 
Signals.  Epicenter  35.  ON  51.  OE  on  February  3,  19>?.  m  =5.0  (C8tGS) 


FIGURE  III-5.  Typical  Half- Amplitude  Rayleigh  Wave  Decay  Patterns,  limes  Indicate  Occurrence 
Of  Half- Amplitude  Increments. 


FIGURE  III-6.  Typical  Half- Amplitude  Rayleigh  Wave  Decay  Patterns.  Times  Indicate  Occurrence 
Of  Half- Amplitude  Increments. 


A/ 16  measurements  were  taken  at  epicentral  distance  less  than  50  degrees. 
Difference  in  signal  dispersion  for  many  paths,  the  wave  period  of  peak  am¬ 
plitude  relative  to  both  source  magnitude  and  period-dependent  propagation 
lactors,  plus  a  variety  of  multipathing  contributions  to  the  main  signal  and 
coda,  all  contribute  to  randomisation  of  the  half-amplitude  times.  Random- 
ntss  in  -he  observations  is  present  in  all  parts  of  the  plots,  however,  with  a 
wide  range  of  decay  times  observed  for  the  A/2  increment,  both  rapid  and 
slow  decay  to  A/4,  and  rapid  and  slow  decay  for  A/8,  etc.  Such  randomness 
might  be  expected  since  we  have  included  a  very  wide  range  of  wave  paths, 
distances,  and  magnitudes  as  a  goal  for  describing  the  generalised  character 
of  the  waveforms  observed  in  a  network  of  recording  stations.  We  could  also 

expect,  and  we  do  see,  definite  non-random  amplitude  decays  for  specific 
source-station  paths. 

Mo  clear  relationship  was  found  for  a  magnitude-half-amplitude 
trend  for  either  network  or  single  path  observations  (magnitudes  range  from 

mb  -  3.2  to  mb  v  6.  2|.  In  fact,  amplitude  decay  patterns  for  specific  paths 

tended  to  be  very  similar  (within  30  to  40  seconds  for  half-amplitude  incre¬ 
ments)  over  at  least  one  order  of  magnitude,  and  paths  repeated  by  signals 
from  about  the  same  magnitude  source  produce  very  nearly  identical  decay 
patterns.  Some  tendency  to  have  longer  first  half- amplitude  (A/2)  decay  time 
with  increasing  magnitude  was  noted  for  specific  paths  (as  was  a  tendency  to 
have  a  longer  wave  period  for  peak  amplitude),  but  this  was  not  a  clear  trend 
in  all  of  the  data  observed.  A  lack  of  very  lar„e  events  (mb  =  7.  0  or  greater) 
prevented  evaluation  of  the  signal  characteristics  of  very  long-term,  multiple- 
path  surface  wave  signals  which  are  well  known  in  observational  seismology 

Distance  distribution  of  the  observations  (number  of  data  channel 

source  pairs)  is  shown  in  Figure  III-8  for  the  entire  data  set.  An  average  of 

about  two  half- amplitudes  were  obtained  from  each  pair,  and  a  maximum  of  four 

times  (A/ 16)  was  the  practical  measurement  limit.  Time  to  noise.  T  .  provided 

n 
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the  possibility  of  a  fifth  measurement,  but  this  is  not  generally  included  in 
the  modeling  unless  the  amplitude  of  noise  equaled  one  of  the  signal  incre¬ 
ments.  Because  of  the  bimodal  characteristic  of  the  data  srmple,  a  sep¬ 
aration  for  analysil  was  made  by  modeling  events  from  either  greater  or  less 
than  50  degrees  of  epicentral  distance. 


Ill-  1  3 


SECTION  IV 


RESULTS 


The  observed  half-amplitude  decay  times  from  the  VLPE  re¬ 
cordings  were  summarized  so  that  the  number  of  observations  with  times 
greater  than  successive  50  second  multiples  are  plotted  and  compared  with 
the  number  of  events  predicted  by  the  gan  ma  distribution  given  in  equation  7. 
The  observed  and  predicted  number  of  events  are  shown  in  Figures  IV- 1,  IV-2, 
and  IV-3.  .'or  the  comparison,  the  gamma  parameter  r  =  0,  and  the  para¬ 
meter  A  =  the  number  of  observations  divided  by  sum  of  all  observed  decay 
times  (or,  equivalently,  the  reciprocal  of  the  mean  decay  time).  The  pre¬ 
dicted  number  of  events  is  the  product  of  the  total  number  of  observations  (N) 
and  P  £\V  £tj,  witb  P  £  W  ^tj  evaluated  at  t  =  50,  100,  150,.  .  . 
seconds.  Differences  in  the  A  term  (or  1/m)  at  each  distance  range  shown 
in  Figures  IV-1,  IV-2,  and  IV-3  are  most  likely  due  to  the  different  data 
bases  used  to  obtain  the  half -amplitude  time  measurements.  Little  duplication 
of  signals  between  the  paper  records  and  digital  seismograms  is  represented 
in  the  observations.  Unprocessed  and  bandpassed  signal  analysis  was  per¬ 
formed  essentially  on  the  same  signal  base,  but  only  60  -  70%  of  the  times 
observed  on  the  bandpassed  recordings  could  be  measured  on  the  unprocessed 
recordings.  Difficulty  in  measuring  second,  third,  and  fourth  half-amplitude 
increments  in  the  wider  band  noise  of  the  unprocessed  recordings  is  the  pri¬ 
mary  cause  of  differences  in  the  latter  two  cases. 


The  data  are  shown  for  signals  observed  at  less  than  or  more 
than  50  degrees  of  epicentral  distance  because  of  the  bimodal  characteristics 
of  the  observations  noted  earlier  (Figure  III-8).  If  the  data  are  shown  to  fit 
a  random  model  in  both  instances,  they  will  also  fit  a  model  which  includes 
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Number  of  Signals  Number  of  Signals 


FIGURE  IV- 1 

Predicted  and  Observed  Gamma  Distributed  Waiting  Times  (Wn) 
For  All  Observed  Half- Amplitude  Rayleigh  Wave  Signal  Amplitude 
Decays.  Gamma  Distribution  Parameters  r  =  0  and  Intensity  1/m. 
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all  observations  if  the  parameter  A  is  changed  appropriately.  Theorems 
regarding  the  randomness  of  subsets  of  observations  taken  from  distributions 
in  the  exponential  class  of  probability  density  functions  (such  as  the  Poisson) 
require  the  foregoing  conclusion.  In  fact,  the  complete  sufficient  statistic 
for  the  entire  data  set  in  each  case  is  At  =  (N^Ajt  +  N2A2t)/  (Nj  +  N2)  lf 
N  is  the  total  number  of  observations  and  the  subscripts  1  and  2  refer  to 
data  observed  at  less  than  or  more  than  50  degrees,  respectively.  Fits  to 
the  total  data  set  for  the  bandpassed  signals  and  unprocessed  signals  (neither 
shown)  for  the  A/ 2  data  alone  using  the  statistic  above  are  satisfactory  at  a 
confidence  level  greater  than  .  950  and  greater  than  .  900,  respectively.  In- 
sufficient  data  to  attempt  a  model  fit  were  found  in  the  paper  record  analysis 
for  events  occurring  beyond  50  degrees  of  epicentral  distance,  which  is  a 
further  reason  for  maintaining  the  separation  of  observed  data  as  shown  here. 

The  "goodness  of  fit"  for  the  models  and  respective  data  sets 
was  tested  by  the  Chi-squared  procedure  in  which  the  predicted  and  observed 
number  of  events  (N)  were  compared  by 

EZ 

(N  pred.  -  N  obs.  )  ^ 

N  pred. 
n 

for  n  of  the  50  second  multiple  intervals.  As  is  customary  for  he  procedure, 

data  for  later  intervals  which  were  rarely  observed  are  lumped  together  in  the 

-  2 

last  (longest  time)  interval.  The  values  of  the  calculated  A  was  then  com¬ 
pared  to  a  standard  table  of  X2  values  for  one  less  degree  of  freedom  than 
the  number  of  intervals  tested.  If  the  calculated  value  of  X  exceeds  the 
tabulated  value,  the  hypothesis  that  the  decay  times  can  be  described  as  a 
Poisson  process  would  be  rejected  at  a  particular  confidence  level.  The  cal¬ 
culations  are  given  in  Tables  IV-1  through  IV- 3,  and  the  results  for  X  are 
given  on  each  of  the  figures. 
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TABLE  IV-1 

GOODNESS  OF  FIT  TESTS  FOR  GAMMA-DISTRIBUTED 
HALF -AMPLITUDE  WAITING  TIMES  (Wn)  IN  RAYLEIGH 
WAVES  RECORDED  ON  UNPROCESSED  DIGITAL  SEISMOGRAMS 


t  (secs) 


N  obs. 

W  >  t 
n 


N  pred. 
W  >  t 


t  (secs) 


N  obs. 
W  >  t 


N  pred. 

W  >t 
n 


53.2 

37.2 


272 

398 

039 


.  182 
.  800 
.495 
.  307 
.039 
.  333 
.  575 


56.9 
47.  6 


.  169 
.406 
.  934 
.  663 
.  158 
.518 
.  321 
.  030 
.496 


Total  Signals  N  =7 

Mean  Waiting  Time  m  =  14 

Distance:  less  than  50  degrees 


greater  than  50  degrees 
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GOODNESS  OF  FIT  TESTS  FOR  GAMMA- DISTRIBUTED 
HALF- AMPLITUDE  WAITING  TIMES  (Wn)  I  I  RAYLEIGH 
WAVES  RECORDED  ON  BANDPASSED  DIGIT AL  SEISMOGRAMS 


w 

n 

N  obs. 

N  pred. 

W 

N  obs. 

N  pred. 

t  ( secs ) 

W  >  t 
n 

W  >  t 
n 

X2 

n 

t  (secs) 

W  >t 
n 

W  >t 
n 

X2 

50 

92 

87.  0 

.288 

50 

87 

85.2 

.  039 

100 

68 

63.  8 

.278 

100 

76 

69.0 

.710 

150 

50 

46.  6 

.248 

150 

65 

57.4 

.  001 

200 

33 

34.  0 

.  029 

200 

51 

47.6 

.242 

250 

24 

25.  0 

.  040 

250 

41 

39.6 

.  050 

300 

17 

18.  1 

.  067 

300 

34 

32.  9 

.  037 

350 

13 

13.  3 

.  006 

350 

28 

27.4 

.  013 

400 

9 

9.7 

.  051 

400 

22 

22.7 

.  022 

450 

6 

7.  1 

.  171 

450 

19 

18.9 

.  000 

>500 

4 

5.2 

.277 

500 

14 

15.7 

.  184 

550 

10 

13.  0 

.  692 

>600 

8 

10.8 

.725 

x2  = 

9 

1,  455 

x2  = 

A11 

3.715 

Total  Signals 

N 

119 

N 

100 

Mean 

Waiting  Time  m  = 

160 

m  = 

270 

Distance:  less  than  50  degrees 

greater  than  50  degree; 

TABLE  IV -3 

GOODNESS  OF  FIT  TESTS  FOR  GAMMA- DISTRIBUTED 
HALF- AMPLITUDE  WAITING  TIMES  (Wn)  IN  RAYLEIGH 
WAVES  RECORDED  ON  PAPER  RECORD  SEISMOGRAMS 


w 

N  obs. 

N  pred. 

n 

t  (secs) 

W  £t 
n 

W  £t 
n 

50 

58 

56.  6 

100 

38 

43.  0 

150 

28 

32.5 

200 

19 

24.  7 

250 

15 

18.7 

300 

13 

14.  1 

350 

11 

10.  7 

400 

9 

8.  1 

450 

8 

6.  1 

500 

7 

4.  7 

550 

5 

3.  5 

600 

3 

2.  7 

650 

2 

2.  1 

>700 

1 

1.5 

13 

Total  Signals 


.  035 
.  582 
.624 
1.310 
.733 
.  086 
.009 
.  010 
.  593 
1.  125 
.  642 
.  034 
.  005 
.  167 

=  5. 955 

=  76 

=  180 


Mean  Waiting  Time  m  = 


Distance:  Less  than  50  degrees 


TT 


Ail  of  fhe  comparisons  of  predicted  and  observed  waiting  time 
distributions  show  a  goodness  of  fit"  value  which  leads  to  a  "do  not  reject 
for  the  hypothesis  that  the  are  Poisson  distributed  at  a  .  90  confidence 
level  or  greater.  Confidence  levels  in  the  results  range  from  about  .925  in 
the  paper  record  analysis  to  greater  than  .  995  in  the  bandpassed  signal  ana¬ 
lysis  for  events  at  less  than  50  degrees.  We  therefore  feel  that  the  data  are 
sufficiently  random  in  a  Poisson  sense  to  proceed  on  with  further  analysis  of 
the  signal  decay  characteristics  in  this  model. 

The  valuable  characteristic  of  observing  the  half-amplitude 
times  as  we  have  and  finding  a  reasonable  fit  to  the  Poisson  model  is  that,  by 
being  able  to  describe  the  half -amplitude  rates,  the  A/4,  A/8,...  rates  should 
also  be  predicted  by  the  gamma  distribution  method.  If  the  time  distribution 
of  the  first  half-amplitude  (A/2)  is  fit  by  the  gamma  parameter  r  =  0,  then 
the  waiting  time  from  peak  amplitude  to  second  half  amplitude  (A/ 4)  should  be 
fit  by  r  =  1,  and  so  on,  without  change  in  the  other  parameters. 

A  demonstration  of  this  characteristic  is  given  for  three  cases 
where  sufficient  times  between  peak  amplitude  A  and  A/ 4  and  A/8  were  ob¬ 
served  to  attempt  a  comparison.  The  data  are  shown  in  Figures  IV-4  and  IV-5 
for  A/4  and  A/ 8  amplitude  increments,  respectively.  Figure  IV-4  shows  the 
observed  distribution  of  the  T^  measurement  for  the  signals  recorded  at  less 
than  50  degrees  epicentral  distance.  The  calculated  number  of  events  r  pected 
was  calculated  from  the  gamma  distribution  using  the  parameter  r  =  1,  and 
X  =  1/m,  with  m  given  in  Figures  IV-1,  IV-2,  and  IV- 3  for  the  respective 

half- amplitude  rates.  The  X  statistic  was  also  computed  for  the  "goodness 
of  fit"  for  each,  and  in  each  case  the  confidence  in  not  rejecting  the  hypothe¬ 
sized  model  is  greater  than  .  950  (see  Tables  IV-4,  IV-5). 

Figure  IV-5  shows  the  observed  and  predicted  distribution  of  the 
number  of  events/waiting  times  for  the  T3  measurement  (time  between  peak 
amplitude  and  one-eighth  amplitude)  for  the  unprocessed  and  bandpassed  signal 
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FIGURE  IV-4 

Predicted  and  Observed  Gamma  Distributed  Waiting  Times  (Wn) 
f  or  One-Quarter  (A/4)  Rayleigh  Wave  Amplitude  Decay  From  Peak 
Signal  Amplitude.  Gamma  Distribution  Parameters  Are  r  =  1  and 
Intensity  1/m.  Epicentrai  Distance  Less  Than  50  Degrees. 
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Predicted 


Observed 

7 


UNPROCESSED  DIGITAL 
N  =  13 

m  =  140  seconds 
*5=  0. 842 

Confidence  level  of  fit  (CF) 
.  950  <CF  <.975 


Wn  >  100  200  300  400  500 

Seconds 


Predicted 


Observed 


BANDPASSED  DIGITAL 
N  =  22 

m  =  160  seconds 
*7=  1.630 

Confidence  level  of  fit  (CF) 
,  975  <CF  <.  990 


Wn£l00  200  300  400  500  600  700  800 

Seconds 

FIGURE  IV -5 

Predicted  and  Observed  Gamma  Distributed  Waiting  Times  (Wn) 
For  One -Eighth  (A/8)  Rayleigh  Wave  Amplitude  Decay  From  Peak 
Signal  Amplitude.  Gamma  Distribution  Parameters  Are  r  =  2  and 
Intensity  1/m.  Epicentral  Distance  Less  Than  50  Degrees. 
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TABLE  IV -4 


GOODNESS  OF  FIT  TESTS  FOR  GAMMA-DISTRIBUTED 
ONE-QUARTER  AMPLITUDE  WAITING  TIMES  (Wn)  IN 
RAYLEIGH  WAVES  RECORDED  AT  LESS  THAN 
50  DEGREES 


w 

n 

t  (secs 

N  obs. 

W  >t 
n  ** 

N  pred. 

W  >t 
n 

X2 

W 

n 

t  (secs) 

N  obs. 

W  >  t 
n 

N  pred. 

W  >t 
n 

X2 

50 

26 

25.  65 

.  005 

50 

40 

39.40 

.  009 

100 

21 

22.  68 

.  124 

100 

33 

35.71 

.205 

150 

15 

19.  14 

.895 

150 

28 

31.49 

.  389 

200 

13 

15.69 

.  461 

200 

24 

26.40 

.  218 

250 

10 

12.  64 

.  551 

250 

20 

22.  06 

.  192 

300 

9 

9.  91 

.  083 

300 

16 

18.  16 

.257 

350 

7 

7.  75 

.  073 

350 

12 

14.63 

.472 

400 

5 

6.  78 

.467 

400 

11 

11.76 

.262 

450 

4 

4.  54 

.  064 

450 

8 

9.39 

.206 

>500 

3 

3.46 

.  061 

500 

6 

2.46 

.286 

550 

5 

5.82 

.115 

>  600 

3 

4.  59 

.  550 

x2,  ■ 

2.  784 

X2 

X11 

=  .842 

Total  Signals  N  = 

27 

N 

=  41 

Mean  Waiting  Time  m  = 

160 

m 

=  140 

Unprocessed  Digital  Seismograms 

Bandpassed  Digital  Seismograms 

Paper  Record  Seismograms 

W 

n 

t  (secs) 

N  obs. 

W  >  t 
n 

N  pred. 

W  >t 
n 

X2 

50 

24 

23.23 

.  026 

100 

20 

21.43 

.  095 

150 

18 

19.  10 

.  063 

200 

14 

16.  75 

.451 

250 

13 

14.  33 

.  123 

300 

9 

12.  10 

.794 

350 

8 

10.  10 

.437 

400 

7 

8.  35 

.218 

450 

6 

6.  89 

.115 

500 

5 

5.  62 

.  068 

550 

4 

4.  85 

.  149 

>600 

3 

3.74 

.  146 

Total  Signals  N 

x2  = 

2.  685 

24 

Mean  1 

Vaiting  Time  m  = 

180 

TABLE  IV -5 

GOODNESS  OF  FIT  TESTS  FOR  GAMMA-DISTRIBUTED 
ONE-EIGHTH  AMPLITUDE  WAITING  TIMES  (Wn)  IN 
RAYLEIGH  WAVES  RECORDED  AT  LESS  THAN 
50  DEGREES 


-  Unprocessed  Digital  Seismograms  EandpIsI7d  Digital  Seismograms 

Wn  ^  °bs.  N  Pred*  W  In  obs.  I  N  pred.  I 

Msecs)  n“  1  Wn“l  X2  t  lo^o\  ^  H  W  >t  v2 


100 

200 

300 

400 

500 

>600 


12.  55 
10.73 
8.  27 
6.71 
3.  68 
2.  50 


.  024 
.  050 
.  195 
.435 
.  038 
.  100 


.  842 
13 


w 

n 

t  (secs) 

N  obs. 

W  >  t 
n  — 

|  N  pre< 

W  > 
n  “ 

100 

20 

21.47 

200 

18 

19.  07 

300 

14 

15.  69 

400 

1  1 

11. 95 

500 

8 

8.  73 

600 

6 

6.  17 

700 
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analysis.  Insufficient  data  from  the  Paper  Records  were  available  for  fitting 
this  increment.  The  fit  in  this  instance  is  not  quite  as  "good"  as  above,  but 
confidence  that  we  should  not  reject  the  model  fit  is  still  very  high  (above 
.  950).  The  fit  is  made  using  the  gamma  distribution  with  r  =  2  ,  and  A  the 
same  as  in  the  r  =  0  and  r  =  1  cases. 
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CONCLUSIONS 

Demonstration  of  randomness  of  Rayleigh  wave  half -amplitude 
decay  times,  within  constraints  of  wide  spatial  distribution  of  sources  and  re¬ 
cording  stations  (and  at  least  a  moderate  range  of  source  magnitudes),  pro¬ 
vides  a  basis  for  pi  edicting  signal  amplitude  -  time  relationships  on  seismo¬ 
grams.  The  two  parameter  gamma  probability  distribution  describes  the  half¬ 
amplitude  occurrence  times  adequately  in  a  statistical  sense.  The  parameter 
A  ,  which  requires  determination  to  use  the  distribution,  can  be  directly  de¬ 
rived  from  a  set  of  observations  of  the  half- amplitude  times  in  recorded  signals. 

Practical  use  of  the  result  given  here  lies  in  estimation  of  sig¬ 
nal  interference  times  expected  in  a  network  of  recording  stations  during  some 
typical  period  of  operating  time.  Using  the  observed  number  of  sources  occur- 
ing  (or  expected)  during  such  a  period,  distribution  of  the  sources  in  terms  of 
magnitude,  and  an  acceptable  amplitude-distance -magnitude  relationship,  the 
peak  signal  amplitudes  can  be  estimated.  Given  these  amplitudes,  the  time 
rate  of  amplitude  decay  by  half -amplitude  increments  can  be  expressed  at  pre¬ 
selected  probability  levels.  A  continuation  of  the  study  will  be  concerned  with 
this  application. 

A  further  continuation  of  the  work  includes  an  evaluation  of  the 
multiple  source  -  single  station  signal  amplitude  decay  times.  If  the  random¬ 
ness  observed  here  is  primarily  a  function  of  multiple  paths  and  varying  source 
magnitude,  randomness  of  the  Poisson  type  should  also  be  found  in  the  single 
station  recordings  with  a  sufficiently  large  base  of  observational  data.  The 
effects  of  the  distance  distribution  of  active  seismic  regions  relative  to  the  sin¬ 
gle  stations  can  then  be  accounted  for,  and  a  more  accurate  estimate  of  poten¬ 
tial  interference  times  by  Rayleigh  waves  in  the  network  will  result. 
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